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^p H ' Abstract 

^ ■ Little information is known about the polarization of gluons 

■ inside a longitudinally polarized proton. We investigate the sensi- 

■ tivity of photoproduction experiments with both beam and target 
longitudinally polarized to the polarization of the gluon distribu- 
tion in the proton. We study the photoproduction of jets and 
heavy quarks and conclude that they are both sensitive to the 
gluon polarization. 
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1. Introduction 



Since the so-called EMC spin crisis has emerged]!], §, much theoretical work 
has been done||. Central to the debate is the size of the gluon polariza- 
tion inside a longitudinally polarized proton. The only information available 
about the gluon polarization is given through higher order corrections to 
spin-dependent structure functions or through its effect on the evolution of 
polarized quark distribution functions. However by considering reactions 
for which the gluon plays a role at leading order, the sensitivity to its po- 
larization will be increased. Many ways to measure the gluon polarization 
in this more direct fashion have been suggested in polarized deep inelastic 
scattering]!, ||, || and polarized hadron-hadron interactions 0. In this paper, 
the sensitivity to the gluon polarization is studied in photoproduction exper- 
iments where both the photon and the proton are longitudinally polarized. 
The production of both jets and heavy quarks is considered. Although these 
processes have been studied before [§, |J, they are reconsidered here in more 
detail and in light of a possible large gluon polarization. 

In Section 2, leading order sets of polarized distribution functions compat- 



ible with the EMC result and the work of Bodwin and Qiu|IJ are presented. 
Different amounts of gluon polarization are considered. As is well known, 
the photon distribution functions play an important role in the photopro- 
duction of jets. The polarized distribution functions of the photon developed 
in Ref. |TT| are used. Their important features are summarized in Section 3. 
In Sections 4 and 5, the results for the photoproduction of jets and heavy 
quarks are discussed, respectively. Finally, in Section 6 the conclusions are 
reviewed. 

2. Distribution functions of the proton 

The relevant quantities are the helicity difference distribution functions: 

Aqi(x p , Q 2 ) = qt(x p , Q 2 ) - q^(x p , Q 2 ), (1) 



where qt(x p , Q 2 )(qi (x p , Q 2 )) is the probability density for a quark of flavor i 
and momentum fraction x p to have a helicity of the same (opposite) sign as 
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the helicity of the proton. The same expression applies for the anti-quarks 
and gluon. Q 2 is the QCD evolution scale. Assuming an SU(3) symmetric 
sea, the following parametrizations are used at Q 2 = Ql = 4 GeV 2 : 



x p Au v (x p , Ql) = N Uv x^(l-x p ) b ^(l + -f Uv x p ) 
Xp/\d v (xp, Qq) = Nd v Xp dv {\ Xp) dv , 
x p As(x p ,Ql) = N s x a p s (l-x p ) b % 

x p Ac(xp,Ql) = 0, (2) 

where 

N Uv = Au v (Ql)/{P{a Uv ,b Uv + 1) + i u J(a Uv + 1, b Uv + 1)), 
N dv = Ad v (Q 2 )/P(a dv ,b dv + l), 

N s = As{Ql)/P(a s ,b s + l), (3) 
/3(x, y) is the Euler beta function, and Aqi(Ql) is the first moment of Aqi(x p , Qq) 

^Qi(Qo)= I A qi (xp,Q 2 )dx p . (4) 

J 



All the quarks are assumed to be massless. Not enough experimental infor- 
mation is available to fit all the parameters, so that some phenomenological 
input has to be used. A strong constraint is provided by the requirement that 
the helicity difference distribution functions must be smaller than the corre- 
sponding unpolarized (helicity-sum) distribution functions. The parameters 
that are used satisfy this requirement with respect to the set D01.1|P"2" 



The small x p behavior is fixed, following Altarelli and Stirling|J, by taking 
a Uv = Udv =0.8 and a s = 0.7. The high x p behavior is assumed to be similar 
to the behavior of the unpolarized distribution functions: b Uv = 3.7, b dv = 4.7, 
and b s = 8. Once an SU(3) symmetric sea is assumed, the normalization of 
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the valence quarks are fixed by the Bjorken sum rule and results from hy- 
peron beta decayS. Following Ref. |, Au v (Q 2 ) = 0.97 and Ad v (Q 2 ) = -0.28 
are used. The sea normalization can be fixed with the EMC result on the 
spin dependent structure function: 



&(xp, Q 2 ) = \H <$(A qi (x p , Q 2 ) + Afcfo, Q 2 )) (5) 



where e; is the charge of qt. Several groups have suggested that the difference 
between the Ellis- Jaffe |TB|] sum rule and the EMC result for the first moment 



of gi(x p , Q ) could be resolved by a larger higher 

order correction due to a large gluon polarization @, [l4j or by a large sea 
contribution 0, [15]] (the sea contribution was set to zero in the Ellis- Jaffe sum 
rule). Several sets of helicity difference 



distribution functions have been developed based on these ideas [0, 16 



However Bodwin and Qiu]l0j showed that a hard gluonic contribution to the 



first moment of gf{x p , Q 2 ) vanishes if an appropriate regularization scheme 
is used. Thus, the only way to resolve the discrepancy between the Ellis 
Jaffe sum rule and the EMC result seems to be to assume a "large" negative 
sea contribution. This point of view is adopted here and As(Q ) = —0.11 
is taken. When the first moment of g\{x P) Q 2 ) is estimated, the data are 
assumed to be Q 2 independent. The EMC result would prefer a somewhat 
larger a s and larger negatively polarized sea when the x p and Q 2 dependence 
are taken into account. However, as mentioned earlier, the unpolarized dis- 
tribution functions give a strong constraint and limit As(Qq) to about the 
value it was assigned. Note that with these values for the parameters, the 
spin carried by the quarks is close to zero and, as was noted in Ref. ||, this 
still needs to be explained. The last parameter ^ Uv can be fit to reproduce 
the x p and Q 2 dependence of gi(x p , Q 2 ) measured by the EMC collaboration 
as well as possible. The result is 7^=2.54. Because of the large sea con- 
tribution, there is no need for a large higher order correction, and a leading 
order analysis is suitable. This is consistent with the fact that later on in this 
paper leading order expressions are used. With this approach there is no ex- 
perimental constraint on the helicity difference gluon distribution functions 



1 Actually only the first moments need to be SU(3) symmetric for this to be true. 
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@. The following parametrization is used: 

x p Ag(x p ,Q 2 ) = N g x° p 6 (l-x p ) 8 , (6) 



where 



N g = Ag(Ql)/P(0.6,1.8). (7) 



The value of the parameters are such that Ag(x p ) < g(x p ) and that Ag(Ql) 
~ 5 is allowed. Three cases are considered: Ag(Ql) = .5, 3.0, and 5.7, 
respectively labeled set 1, 2, and 3. In Fig.l both the helicity difference and 
unpolarized quark and gluon distribution functions are shown at = 4 GeV 
for the three sets. The distribution can be obtained at another Q 2 by evolving 



them with spin-dependant Altarelli-Parisi equations |17|- In Fig. 2 the EMC 
results, with statistical and systematic uncertainties added in quadrature are 
compared to the results of set 1 and 3, calculated at the appropriate Q 2 . 
The agreement between the data and the sets is good considering the fact 
that the sets are not the result of a fit of all the parameters. The difference 
between the two sets is due to the difference in the evolution of the sea quark 
helicity difference distribution functions. 

In this section, three sets of helicity difference distribution functions were 
developed. At Q 2 = Ql, they differ only in the size of the gluon polariza- 
tion. As was pointed out, even the quark distribution functions are not well 
constrained yet. For the purpose of this paper, it is assumed that the quark 
distribution functions are fixed and the emphasis is on the sensitivity to the 
gluon helicity difference distribution function. 



3. Distribution functions of the photon 



As in the unpolarized case[lS], an "asymptotic" solution can be derived from 



the spin- dependent Altarelli-Parisi evolution equations in the limit of large 



2 The constraint from evolution is small. The situation is similar to a leading order 
analysis of F% in the unpolarized case. 
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momentum fraction (x 7 ) and Q 2 . This is due to the direct coupling of the 
photon to quarks. In that asymptotic limit, the Q 2 dependence of the helicity 
difference distribution functions can be factorized: 



Aq t (x y ,Q 2 ) = ^ln(j^)Ahi(x~t), 



Ag(x„Q 2 ) = ^ln&Ah a (x,), (8) 



where a is the QED coupling constant and A the QCD constant. The Ah 
functions depend only on x 1 . The expression in Eq. |9] can be introduced in 
the spin-dependent Altarelli-Parisi evolution equations and the remaining x 1 
dependent equations for the Ah can be solved numerically. This was carried 



out in Ref. |TT] where a parametrization of the solution is provided 0. Four 
flavors and massless quarks were assumed. A solution valid at lower Q 2 and 
x 7 would require input distributions functions at a fixed Q 2 . Unfortunately, 
there are no experimental data to fix the parameters of these distribution 
functions, so the asymptotic solution is utilized for the entire Q 2 range under 
consideration. The same approach can be used in the unpolarized case and 
the parametrization of Ref. |I] is used here. In Fig. 3 both the polarized 
and unpolarized distribution functions are shown for up-type and down-type 
quarks. The distribution functions for the quarks are harder than in the 
proton case, as can be seen by comparing Fig. 1 and 3. This will favor 
configurations with high x 7 and low x p . The polarized distribution function 
of the gluon is also shown in Fig. 3 even though at the energy considered 
in this paper, it doesn't play any role. Notice that at high x-y the helicity 
difference distribution functions of the quarks are positive. 



4. Two— jet production 

Our goal is to study the sensitivity of the photoproduction of two jets to 
the gluon polarization. Both the photon and the proton are longitudinally 

3 At high x 7 the parametrization used for the quarks is bigger than the numerical result 
and the unpolarized distribution functions. Here, to correct for this, Aq(x 7 ) = q(x 7 ) is 
used for x~, > 0.95. 
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polarized. The photoproduction of two jets receive contributions from two 
classes of subprocesses. In the first class of subprocesses the photon inter- 
acts directly with the constituents of the proton (the "direct" contribution). 
In the second class, the photon interacts through its distribution functions 
(the "resolved" contribution). At leading order in perturbative QCD, the 
integrated cross section for the photoproduction of two jets is given by the 
following expression: 



a, 



7p— >2j'ets I 



^ ; J f a /ry(Xry, Q ) / dX p f b/p (X p ,Q 

a,b < 



&ab^2jets{S) 



(9) 



where s = 2E 1 M p is the square of the total energy available and s = x^x p s 
is the square of the center of mass energy of the subprocesses. f a /~/(fb/p) 
is the distribution function of parton a(b) in the photon(proton). In the 
direct contribution case, /a/7 (^7) is replaced by 5(1 — x-y). The choice of the 
scale Q 2 bears the usual ambiguity, Q 2 = is used. The sum is over 

all the possible subprocesses. The direct contribution is composed of two 
subprocesses: jq — > gq and — > qq. The resolved contribution is composed 
of eight subprocesses. The dominant subprocesses at the energy considered 
in this paper are qq' — > qq' and gq — > gq. The lower limits of integration are 
given by the following expressions: 



T 



S 



T ~< = — ^ 10 ) 



where pTmm is the minimum px of the jets. The matrix elements necessary 
to calculate the subprocess cross sections, a, can be found in Ref. [2(] and 
21 . The integrated cross sections are presented in Table [j], at E 7 — 200 GeV 



and 400 GeV, and PTmin = 3 and 5 GeV. E 1 = 200 GeV corresponds to 
the average value for E 1 of present unpolarized experiments and 400 GeV is 
about the upper limit. pTmin = 3 GeV is the lowest value at which jets have 
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been observed in fixed target experiments p2[] . pTmin = 5 GeV is presented 
to show the variation of the different contributions with PTmin • Both the 
direct and the resolved contributions are presented. The two contributions 
are furthermore divided into quark and gluon contributions corresponding to 
subprocesses involving a quark or a gluon inside the proton, respectively. The 



unpolarized case has already been study extensively [p0| , |23|| . The important 



points are summarized briefly. As can be seen from Eq. [11], both the x p and 
x 7 thresholds decrease when the energy increases or the prmin decreases. As 
a result, the relative size of the gluon contribution (in both the direct and 
resolved contribution) and of the resolved contribution increases whenever 
the energy is increased or the minimum px decreased. 

The integrated helicity difference cross section for the production of two 
jets is given byH: 



(°" ~ ° )jp^2jets 

a,b 



s ) = Y. f dx i A /a/ 7 (a; 7 ,Q 2 ) 

a,b 

( dx p Af b/p (x p ,Q 2 )Aa ab ^ 2jets (s). (11) 



The first sign superscript in cr ++ and cr + ~ corresponds to the helicity of the 
photon, and the second to the helicity of the proton. The A/ a / 7 and Af b / P 
are the helicity difference distribution functions as defined in section 2 and 
3. Acr is the helicity difference cross section of the subprocess. The matrix 



elements for the different subprocesses can be found in Ref. [5] and [IT]. The 
longitudinal asymmetry can now be defined: 



As usual, it is advantageous to consider a ratio because the theoretical and 
experimental uncertainties tend to cancel out. The results for the asymmetry 
of each of the contributions are shown in Table |2| for set 1 and 3 (smallest and 
largest gluon helicity difference distribution function), at the same energies 
and pTmin as in Table |]. An important property of the asymmetry is that 
it is a weighted average of the asymmetries of each of the contributions. 



8 



First, the direct contribution is considered. The quark contribution is 
dominated by the u v quark because the cross section is proportional to the 
square of the charge of the quark, and the Au v is the largest of the quark 
helicity difference distribution functions. Both Au v and Aa iq ^ gq are positive 
so that the quark contribution gives a positive asymmetry. The difference 
between set 1 and set 3 for the quark contribution is small; it is due to the 
difference in evolution of the quark helicity difference distribution functions. 
For the gluon contribution Ag is positive and Aa ig ^ qq is negative such that 
the asymmetry of this contribution is negative. As expected, there is a large 
difference between the results of the two sets in this case. At E 1 = 200 GeV 
and pTmin = 5 GeV the difference between the asymmetry of the two sets 
for the gluon contribution is about 60%, whereas the difference for the total 
asymmetry is only about 15%. The large difference in the gluon contribution 
doesn't survive, because the size of the gluon contribution to the cross section 
is relatively small, see Table [1|. For the other three cases in Table 0, the size 
of the gluon contribution is bigger and the difference between the two sets 
for the total asymmetry of the direct contribution is about 35 to 50%. Notice 
that it is important to consider the difference between the two sets, and not 
just the result of each set separately. 

Second, the resolved contribution is studied. In this case, both the quark 
and gluon contributions give a positive asymmetry. As in the direct case, 
the difference between the two sets for the quark contribution is small. The 
difference between the two sets for the gluon contribution is not as big as in 
the direct case. This is due to the fact that the asymmetry of the leading 
subprocess (gq — > gq) is not as big, and that the helicity difference distribu- 
tion functions of the photon had to be folded in. As a result, the difference 
between the two sets for the total asymmetry of the resolved contribution is 
at most 15%. 

As can be seen in Table [|, the difference between the two sets for the 
total asymmetry is of the order of 10-15%. This difference is rather small 
compared to the range span by the gluon contribution of the direct con- 
tribution. The problem stems from the fact that the gluon contribution is 
negative in the direct case and positive in the resolved case, such that the 
two contributions partially cancel each other. An obvious way to improve 
upon this is to separate the direct and resolved contributions, and then use 
the direct contribution to measure the gluon polarization, as it is the most 
sensitive contribution. Eventually, the resolved contribution could be used 



9 



to study the quark helicity difference distribution functions of the photon. 
The same techniques developed for the unpolarized case can be implemented 
to separate the direct and resolved contributions. One way to do this is by 
tagging the remnant jet coming from the photon in the resolved case [p4 



an- 



other is by a complete reconstruction of the kinematics and an appropriate 
cut on a; 7 [^5] . 



More detailed information can be obtained by looking at the longitudinal 
asymmetries of the differential cross section. In Fig. 4 the ^-distributions 
are presented at E 1 = 200 GeV and pTmm = 3 GeV for the three sets de- 
veloped in Section 2. In Fig. 4a the unpolarized cross section is presented 
for the direct (dashes) resolved (dots), and total (solid) contributions. The 
resolved contribution peaks at higher x p to compensate for the lower x 1 . In 
Fig. 4b, 4c and 4d the asymmetry for the direct, resolved and total contribu- 
tions, respectively, are plotted for the different sets (dashes: set 1, dots: set 
2, solid: set 3). The weighting between the direct and resolved contribution 
to form the total contribution is well apparent. For each contribution, the 
largest difference between the three sets is at low x p . Note that it is possible 
to reconstruct x p only if both jets in the event are measured. In case this can 
not be done, the rapidity distribution of the jets in the 7-proton center of 
mass is shown in Fig. 5. Positive rapidity is in the direction of the incoming 
photon. The largest difference is at larger rapidity that corresponds to the 
lowest x p . One could also look at the pt distributions, and the difference 
between the sets is similar. For comparison, in Fig. 6 and Fig. 7 the same 
plots as in Fig. 4 and Fig. 5, respectively, are shown for E 1 = 400 GeV and 
PTmin = 3GeV. 



5. Heavy Quark production 

The formulas presented in Section 4 for the photoproduction of two jets 
are also valid for the photoproduction of an heavy quark pair. The only 
modification is the replacement of PTmm by m q , the mass of the heavy quark, 
in Eq. 10. We will consider the production of the charm quark with m c = 
1.5 GeV. As is well known, the resolved contribution for the photoproduction 
of heavy quarks for the energy range considered here is of the order of a few 
percents, and can be neglected. For the direct contribution there is only one 
subprocess at leading order in perturbative QCD: jg — > QQ, where Q stands 



10 



for a heavy quark. Assuming that there is no spin effect in the fragmentation 
of the charm into a D-meson, the Peterson fragmentation function can be 
used pfifl : 



D^{z) = Nz{l-zf/{{l-zf + ezf (13) 

where z is the momentum fraction of the D-meson, and N is taken such that 
D®(z) is normalized to 1. mo ~ Tn c is assumed as in the derivation of Eq. 14. 
The parameter e is taken at .15. The matrix elements needed to calculate 
the asymmetry were evaluated using the method described in Ref. f?7|. The 
results for the integrated cross section and the asymmetry are presented 
in Table || for sets 1 and 3. The difference between the asymmetries of 
the two sets is about 12% at E 1 = 200 GeV and 5% at E 1 = 400 GeV. 
These results are actually misleading. In Fig. 8 the p^-distribution of the 
D-meson along with the asymmetry distribution are presented for E 7 = 
200 GeV. It is apparent that at low px where the mass terms dominate in 
the cross section, the asymmetry is positive, whereas at high px where the 
mass terms are not as important, the asymmetry is negative (as in the two jet 
production case). The difference between the sets is bigger than suggested 
by the integrated asymmetry. A better variable in this case to describe the 
difference between the sets is given by an "absolute" asymmetry, \A\u, where 
instead of integrating the asymmetry at each phase space point with its sign, 
the absolute value at each point is integrated. The difference between the 
two sets for the absolute asymmetry give a measure of the biggest difference 
that can be reached. The results for \A\u are also shown in Table [| The 
difference between the sets is of the order of 20% at E 1 = 200 GeV and 15% 
at Ery = 400GeV. The lower energy is favored in this case. In Fig. 9 and 
Fig. 10 the rapidity distribution and sum of the rapidity distributions (in 
case both D-mesons are reconstructed) are presented. As it is unlikely that 
the kinematics of the whole event can be reconstructed, the x p distribution 
is not shown. 



11 



6. Conclusions 



We have shown that both jets and heavy quark production in a photoproduc- 
tion experiment can lead to a successful measurement of the gluon helicity 
difference distribution function. Considering the total asymmetry, the two 
jets and heavy quarks production cases have similar sensitivity. However, 
the small difference in the asymmetries between the different gluon polar- 
izations (10-20%) might be a limiting factor. The best way to measure the 
gluon helicity difference distribution function is by using two jet production 
at low px , with separation of direct and resolved contribution. The direct 
contribution has the biggest sensitivity, with differences in the asymmetries 
of the order of 35-50%. The resolved contribution could be used to study 
the quark helicity difference distribution functions of the photon. 

To conclude we emphasize that a direct measurement of the gluon helicity 
difference distribution function is important as it would clarify some of the 
theoretical debate. 
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Tables 







direct 


resolved 


tot 




PT min 


q-cont 


g-cont 


tot 


q-cont 


g-cont 


tot 




(GeV) 


(GeV) 


(nb) 


(nb) 


(nb) 


(nb) 


(nb) 


(nb) 


(nb) 


200 


5 


10.2 


3.2 


13.4 


2.5 


1.2 


3.7 


17.1 


200 


3 


103. 


110. 


213. 


91. 


110. 


201. 


414. 


400 


5 


23.7 


17.3 


41.0 


13.9 


12.6 


26.5 


67.5 


400 


3 


129. 


217. 


346. 


231. 


371. 


602. 


973. 



Table 1: Integrated cross sections for dijet production averaged over the 
initial spins for different E 7 , and pTmm • 









direct 


resolved 


tot 


set 




PT min 


q-cont 


g-cont 


tot 


q-cont 


g-cont 


tot 






(GeV) 


(GeV) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


1 


200 


5 


34.8 


-7.9 


24.7 


10.1 


3.4 


7.9 


21.1 


3 


200 


5 


34.9 


-74.5 


9.1 


10.1 


30.1 


16.8 


10.7 


1 


200 


3 


24.9 


-8.3 


7.8 


2.8 


2.6 


2.7 


5.3 


3 


200 


3 


24.9 


-93.2 


-36.0 


2.8 


29. 


17.2 


-10.1 


1 


400 


5 


28.4 


-8.9 


12.7 


4.7 


3.2 


4.0 


9.3 


3 


400 


5 


28.9 


-89.2 


-20.8 


4.8 


31.6 


17.5 


-5.7 


1 


400 


3 


18.0 


-7.5 


2.0 


.9 


1.8 


1.4 


1.6 


3 


400 


3 


18.2 


-84.6 


-46.4 


.9 


20.3 


12.9 


-8.7 



Table 2: Asymmetries for dijet production for different sets of polarized 
distribution functions (set 1 and 3 of section 2), E 7 , and prmin • 
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proton 




cross 


An 


\A\u 


set 




section 








(GeV) 


(nb) 


(%) 


(%) 


1 


200. 


993. 


1.2 


2.2 


3 


200. 


993. 


13.2 


23.6 


1 


400. 


1475. 


.5 


1.5 


3 


400. 


1475. 


5.8 


16.7 



Table 3: Integrated cross sections and asymmetries for heavy quarks pro- 
duction for different sets of polarized distribution functions (set 1 and 3 of 
section 2), and E T 
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Figure 1^ 
difference "~ 
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^tr^^^n^i^g fesji[^) n ^§ th^r^oji Q 2 Q = AGeV 2 : a) 

k, c) sea quark, and d) gluon. In Fig. 
d) the 3 dashi^d ckrves correspond frori the lowest to the highest to set 1, 2, 
and 3. 



Figure 2: oE^fC- result -for-trle^^pin cependant structure function g\ as a 



function of St ^-Mso shown a^c ' thcn^CCT ilts for set 1 (dashes) and set 3 (solid 
line), calculated at the appropriate Q 2 . 



Figure 3: Unpolarized distribution functions (Ref. 0, solid) and helicity dif- 
ference distribution functions (Ref. ITTL dashes) of the photon at Ql = 4GeV 2 : 
a) up type quark, b) down type quark, c) gluon. 
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Figure 10: Same as in Fig. |8] for the rapidity-sum distribution. 
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